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Abstract 

The recent observation of the purely leptonic decay — ► t^v t at the B factories per- 
mits a sizeable contribution from a charged Higgs boson (H^). Such a would also 
contribute to the decays Df — ► and — ► t^v t , which are being measured with 

increasing precision at CLEO-c. We show that the branching ratios of Df — > fi v a and 
Df — ► t^v t could be suppressed by up to 10% from the Standard Model prediction, 
which is larger than the anticipated precision in the measurements of these decays at 
forthcoming BES-III. 
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1 Introduction 



In April 2006 the BELLE collaboration announced the first observation of the purely leptonic 
decay 5 ± — > t ± v t [1] utilizing an integrated luminosity of 414 fb _1 . The measured branch- 
ing ratio (BR) is in agreement with the Standard Model (SM) rate within theoretical and 
experimental errors: 

B R(£?± - T ±v T ) = (1.79t° it(stat)t° fi(syst)) x 10" 4 (1) 

Subsequently, BABAR reported an improved upper limit with 288 fb _1 [2]: 

BR(B ± -> t±v t ) < 1.8 x 10~ 4 90% ci (2) 

Significantly improved precision for BR(_B ± — > t^v^) would require a high luminosity upgrade 
of the existing B factories [3]. Such a facility would also provide sensitivity to the SM rate for 
the analogous leptonic decay — > /i ± z/ /J . 

In the SM the decay 5 ± — > t^v t proceeds via annihilation of to a W*. It has been 
known for some time that New Physics could significantly affect BR(5 ± — > r ± z/ T ). A charged 
Higgs boson (i^ ± ) present in the Two Higgs Doublet Model (2HDM) would also mediate the 
annihilation of B ± [4], and for Model II type Yukawa couplings its effect on BR(_B ± — > r ± z/ T ) is 
determined by the ratio R = t&nj3/m H (where mu is the mass of H ± , tan/3 = v 2 /vi and Vi,v 2 
are vacuum expectation values of the Higgs doublets). The contribution from interferes 
destructively with the W ± mediated SM diagram, and thus the experimental measurement of 
BR(£? ± — > t^Vt) leads to two allowable regions for R [5]: 

(i) A region where the contribution is small (a perturbation to the SM rate), roughly 
corresponding to < R < 0.15 GeV -1 . 

(ii) A region where if ± contributes sizeably but the large destructive interference maintains 
the SM rate, roughly corresponding to 0.22 GeV" 1 < R < 0.33 GeV -1 . 

In a more general 2HDM in which the bottom quark receives mass from both Higgs doublets 
(e.g. which includes the Minimal Supersymmetric SM at the one-loop level) the above regions 
for R are altered [6]. Of particular interest is region (ii) since such sizeable values of R might 
have an observable effect on other processes in which H ± contributes at tree-level. 

The hitherto unobserved exclusive semi-leptonic decay B — > (where B can be charged 

or neutral) has a relatively large BR (ps 8 x 10~ 3 ) in the SM, and can also be mediated 
by H ± [7], [8]. A first observation of this decay would require several ab _1 of integrated 
luminosity. The angular distribution of the r [9] would also provide information on R with 
reduced uncertainties from Cabbibo-Kobayashi-Maskawa (CKM) matrix elements and strong 
interaction effects. Direct production of H ± at the Tevatron (via t — > i/ ± 6) [10] and at the 
LHC (e.g. via gb — > tH^) would also probe the region (ii). Notably, LHC simulations [11] 
promise verification or falsification of region (ii) for tan/3 < 100, as well as probing a sizeable 
fraction of region (i). 

In this paper we consider the effect of H ± on the purely leptonic decays of the charmed 
meson Df. Due to the smaller Yukawa couplings of H ± to the lighter quarks the effect of H ± 



on such decays will clearly be smaller than the effect on B^ 1 — > t^Vt. However, BR(Df — > 
V^V/jt) ~ 5 x 1CT 3 and BR(D^ — > t^Vt) ~ 5 x 10~ 2 are much larger than BR(B ± — > r ± i/ r ) 
and thus small perturbations to the SM rate from might be accessible to ongoing and 
forthcoming experiments. Measurements of BR(D f -> /i ± z/ M ) and BR(Df — > r ± z/ T ) have been 
performed by a variety of experiments (for a summary see [12]) and are being improved by 
ongoing CLEO-c [13], [14]. Moreover, BABAR has recently observed Df — > y^v^ [15], with 
a precision comparable to that at CLEO-c. At BES-III, which is due to commence in late 
2007, a precision of 2% for BR(D^ — > /i ± ^) and BR(D^ — > t^t) is expected after 4 years of 
operation at design luminosity [16]. 

To date, the main motivation for searching for Df — > /i ± z/^ and Df — > r^Vr has been to 
measure the decay constant fo s , assuming that no New Physics contributes. Such measurements 
are vital tests of lattice predictions for the heavy quark systems. However, if if ± is contributing 
sizeably to B — > r^Ur it would have a non-negligible effect on Df — > \i v ^ and Df — ■> r v T . 
This was noted in [4], [17] and a first quantitative study was performed in [18]. Motivated 
by the recent observation of B — > r ± z/ r , in this paper we update and extend the analysis 
of [18], assuming that is contributing sizeably to B — > r v T . We point out that a 
precise measurement of the ratio BR(_D^ — > yU ± z/ At )/BR(D ± — > /i ± z/ /J ) (which will be possible 
at BES-III) could provide a hint for the presence of H ^ . Our work is organized as follows. In 
section 2 we review the contribution of if ± to the leptonic decays and summarize the current 
experimental situation. The numerical analysis is contained in section 3 with conclusions given 
in section 4. 



2 The decays B ± / ± z/ / and Dt, / ± z/ / 

In the SM the purely leptonic decays B ± — > l^vi and — > l^vi proceed via annihilation of 
the heavy meson into W*. Singly charged Higgs bosons, which arise in any extension of the 
SM with at least two SU(2)l x U(1)y Higgs doublets with hypercharge Y — 1, would also 
contribute to these decays. For B~ — > l~V\ the Feynman diagram is as follows: 



B- 



The tree-level partial width is given by [4]: 
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where Gp is the Fermi constant, mi is the mass of the lepton, mg is the mass of the B meson, 
V u b is a CKM matrix element, and fg is the decay constant. In the 2HDM (Model II), in which 
the b quark only couples to one of Higgs doublets at tree-level, the scaling factor r# of the SM 
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rate is given by: 

r H = [l-ml t -^l}^[l-mlR"f. (4) 

The H ± contribution interferes destructively with that of W ± . There are two solutions for 
r H = 1 which occur at R — and R = 0.27 GeV -1 (see Fig{Tk). If the b quark couples to 
both Higgs doublets at tree-level, which is referred to as the 2HDM (Model III) [19], Eq. (j3J) is 
modified to [6]: 

/ tan 2 /5 m 2 B \ 2 

r H=[l- —2~ ) ( 5 ) 



1 + e tan (3 rn 2 H± 

In the Minimal Supersymmetric SM (MSSM) the parameter Iq is generated at the 1-loop level 
[20], [21] (with the main contribution originating from gluino diagrams) and may reach values 
of 0.01. The redefinition of both the b quark Yukawa coupling and the CKM matrix element 
Vub are encoded in e [22]. The impact of e ^ on r H has been developed in [8], [9], [23], [24]. 
In particular, the value of R where r# = and r# = 1 shifts depending on the magnitude and 
sign of io- 

For the decay — > l~V\ the partial width is given by [4]: 

m = g|mC l"" 2/ °'-' |^,r (l - x r w (6) 



57T \ m 



where mr, (g) is the mass of the meson, V C( i( cs ) are CKM matrix elements, fn (s) are decay 
constants, and rv s ) is the analogy of r# given by: 

r w = [1 - ml ^l { ^^ )] 2 ^ {1 _ m 2 /2 2 ( — — — )] 2 (7) 

As in the case for B ± — > / ± ^, the H ± contribution depends on R = tan /3 /mjf±, although 
the factor m 2 D m q j[m c + m q ) in Eq ,([7j) is considerably smaller than mg present in Eq .(T4|). 
The H ± contribution to — > l^vi is essentially negligible due to the smallness of m^/m c and 
thus r ~ 1. However, for Df the scaling factor r s may differ from 1 due to the non-negligible 
fn s / m c [4], [17], [18]. The expression for r s in Eq. ([7]) can contain an eo correction analogous 
to Eq. (j3J), but we do not write this contribution explicitly due to the considerable uncertainty 
in m s j (m c + m s ). In Fig{T]D we plot r s as a function of R. Although the effect of H ± is only a 
perturbation to the SM rate for BR(Df — > h ± i i ^ 1 t ± i> t ) ) such deviations might be measurable 
since their BRs in the SM are much larger than BR(B ± — > t^v^. Prospects for precise 
measurements of BR(D^ — > fi^Vfj,, f^v^) are bright. The ongoing CLEO-c programme [27], [28] 
(which will terminate around April 2008) utilizes the production mechanism e + e~ — > Df*Df, 
and expects a final precision ~ 10% for BR(Df — > /j, ± i> fl ,T ± i' T ). At forthcoming BES-III an 
integrated luminosity of 20 fb _1 (which corresponds to 4 years at design luminosity) would 
enable a precision ~ 2% for BR(D^ — > p^v^r v) [16]. In Table 1 we show the SM rates 
and the current measurements of the purely leptonic decays D ± — > I v\ and Df — > l^vi from 
CLEO-c. The displayed measured value of BR(Df — > r^v) is an average of separate searches 
for — > and — > e ± vv. For our values for the SM BRs we take fjj = 200 MeV [29] 
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Figure 1: Scaling factor of the SM rate (solid line) as a function of R(= tan/3/m^±). Left 
panel (Fig. la): r H for B ± — > t 1 * 1 //,-; Right panel (Fig. lb): r s for £>± — > h ± v^ 1 t ± v t with 
m s / (m c + m s ) = 0.08. 



Decay 


SM BR 


Current Exp BR 


CLEO-c error 


BES-III error 


D ± -> e ± z/ 


8.24 x 10~ 9 


< 2.4 x 10~ 5 [25] 








3.50 x 10" 4 


4.40 ± 0.66t£™ x 10~ 4 [25] 


~ 10% 


-2% 




9.25 x 10~ 4 


< 2.1 x 10~ 3 [26] 






L>± -> e ± z/ 


1.23 x 10~ 7 


< 3.1 x 10~ 4 [13] 






£>± -> n±u 


5.22 x 10~ 3 


6.57 ±0.9 ±0.28 x 10~ 3 [13] 


~ 10% 


-2% 


Df ^ r ± z/ 


5.09 x 10~ 2 


6.5 ±0.8 [14] 


~ 10% 


~ 1.5% 



Table 1: The leptonic decays of D ± and Df: BR predictions in the SM (uncertainty ~ ±30% 
not shown), current measurement from CLEO-c (third column), expected final precision at 
CLEO-c (fourth column) and BES-III (fifth column). 

and fo s = 250 MeV [29], [30]. The current uncertainty of ~ 15% in the lattice calculations 
of the decay constants induces an error of ~ 30% for all leptonic BRs. Previous searches for 
D ± — > fx^u^ can be found in [31] and a summary of measurements for BR(D^ — > i^v^r ± v 1 ) 
prior to the CLEO-c programme is given in [12]. 

In addition to the above charm facilities, the large amount of charm data at the B factories 
can be used to search for leptonic decays of D ± and Df. At BABAR, a search for the process 
Df* — > Df^f — > /i ± z/ M 7 was performed with 230 fb~ 4 of integrated luminosity, with measured 
BR [15]: 

BR(Df -> fj^u) = (6.74 ± 0.83 ± 0.26 ± 0.66) x 10~ 3 (8) 

The third error (which is not present in the CLEO-c measurement) is from BK(Df — > 07r ± ), but 
should decrease with improved measurements of BR(_D^ — > (pn^) at CLEO-c and forthcoming 
BES-III. No search for Df — > /i ± z/^ has yet been performed by BELLE. In summary, the 
precision for the measurements of BR(Df — > h ± i , ^ 1 t ± u t ) is constantly improving and will 
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Figure 2: Contours of r s in the plane [R,m sc ], where R = tan/3/m#± and m sc = m s j '(m s + m c ). 

approach the percent level in the next few years. In contrast, a significant improvement in 
the measurement of BR(5 ± — > t^Vj) would require a luminosity upgrade of the existing B 
factories. 

3 Numerical Results 

The current measurement of — > r ± z/ T (Eq. p])) restricts i? to an interval which can be used 
to predict a region for r s i.e. the deviation of BR(Z) ± — > // i^, r ± i/ T ) from the SM rate. The 
BELLE collaboration [1] uses V u b = 4.39 ± 0.33 x 1CT 3 and the unquenched lattice result = 
216 ± 22 MeV [33] which give the theoretical prediction BR( J B ± -> r ± z/ T ) = 1.59 ± 0.4 x 10~ 4 . 
Comparing this value for BR(S ± — > t u T ) with the measurement in Eq. ([1]) gives the scale 
factor r# = 1.13 ± 0.51. The interpretation in the 2HDM leads to two allowable regions for R. 
The region of larger R is 0.22 GeV" 1 < R < 0.33 GeV" 1 for e = in Eq. p. Note that using 
smaller values of and V u b for the SM prediction for BR(5 ± — > r 1 * 1 //,-) would shift the above 
allowed region for r# to larger values. Therefore for our evaluation of r s we consider 0.20 GeV -1 
< R < 0.40 GeV -1 . For the quark masses m s and m c we use the Particle Data Group [32] values 
and obtain 0.03 < m s /(m c + m s ) < 0.12 at the 2a level. However, we will consider the range 
0.03 < m s /(m c + m s ) < 0.15 which allows for enhancement of the s quark Yukawa coupling 
(which is the source of the numerator in m s /(m c + m s )) in the context of the 2HDM (Model 
III). In Fig. [2]we plot contours of r s in the plane [R, m s j (m c + m s )]. The suppression of the SM 
rate is between 2% and 15% in most of the parameter space. Thus the contribution can be 
larger than the expected BES-III precision (~ 2%) for BR(D^ — > /i ± ^ M , r ± z/ T ) and comparable 
to the anticipated error in future lattice calculations of fn a - The presence of H ± would lead to 
a deceptively smaller measured value of the decay constant f Dg . We suggest that the possible 
effects of any should not be overlooked when comparing the experimentally extracted fr> a 
to the lattice calculations. Importantly, discovery of in direct production channels at the 



\ 






\o.85 












\0.95 






: V =0-9 


8 















5 



16 




R 



Figure 3: TZ^ as a function of R(= tan f3/m H ±) for m s /{m s + m c ) = 0.08 and f D J f D =1.22, 
1.24 and 1.26. 

Tevatron or LHC with properties corresponding to the region of 0.20 GeV -1 < R < 0.40 GeV -1 
would necessitate inclusion of the scaling factor r s in BR(.D^ — > fi^is, t^v^) in order to do a 
proper comparison of lattice calculations with experiment. 

The sizeable error in the decay constant makes it difficult to confirm the presence of 
from measurements of BR(-D^ — > y^v, r u T ) alone. An additional observable which is more 
sensitive to the contribution from H ± is the ratio of the muonic decay rates IZ^ defined by 

K, = BR(Df - f i ± u,)/BR(D ± - ^) ~ r s (f D Jf D ) 2 (9) 

The effect of H on BR(Z) =t — > // u^) is negligible and TZ^ is proportional to r s . The lattice 
prediction for Jds/Jd is known with greater precision than the prediction for the individual 
values of the decay constants and hence TZ^ is a better probe of the H ± effect on BR(.D^ — > 
fj^u^). A similar ratio (7Z T ) for the tauonic decays could be measured at BES-III [16] (and 
possibly at CLEO-c) although we only consider TZ^ for which both experiments expect superior 
precision. An unquenched lattice calculation [29] gives JdJId = 1.24 ± 0.07 ± 0.01. The 
first measurements of f D Jf D from CLEO-c [14] and BABAR [15] are 1.26 ± 0.11 ± 0.03 and 
1.27 ±0.14, where both use the CLEO-c result for fo [25]. The central values are in agreement 
with the above unquenched lattice result, and differ from the quenched results which typically 
give f D Jf D = 1.10 [34]. 

In Fig. [3] the ratio TZ^ is plotted as a function of R for m s /{m s + m c ) = 0.08. We take 
fDslfn — 1.24 ±0.02 which might correspond to the error in the lattice calculations as BES-III 
approaches 2% — 3% precision for these decays. For fnj Id = 1-24 one can see that 7Z^ varies 
from 14.7 to 13.3 (i.e. a 10% suppression) as R varies from — > 0.4 GeV -1 . Hence precise 
measurements of TZ^ could favour or disfavour H^. 
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4 Conclusions 



The recent observation of B — > does not preclude a sizeable contribution from H in 
the 2HDM or MSSM. Such a H would contribute to the leptonic decays Df — > fi u^, r u T . 
We showed that H ± can suppress BR(D^ — > h ± i> im1 t ± i> t ) by 10% or more, which is larger 
than the expected experimental error (2% — > 3%) for future measurements of these decays. 
We emphasize that new physics effects like these should not be overlooked when compar- 
ing the experimental measurements of fu s to the lattice predictions. Moreover, favouring 
or disfavouring the presence of could be obtained by precise measurements of the ratio 
BR(Df — > / u ± z/ At )/BR(D ± — > fJ^v^). Prospects for Df — > p^v^ and Df — > f^v T are bright 
with ongoing CLEO-c and forthcoming BES-III, and we encourage their study together with 
B ± — > r ± z/ at the B factories. 
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